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Abstract 

^ I Frictional sliding along an interface between two identical isotropic elastic plates 

0^ ' under impact shear loading is investigated experimentally and numerically. The 

plates are held together by a compressive stress and one plate is subject to edge im- 
\^ ] pact near the interface. The experiments exhibit both a crack-like and a pulse-like 

CD I mode of sliding. Plane stress finite element calculations modeling the experimental 

configuration are carried out, with the interface characterized by a rate and state 
dependent frictional law. For low values of the initial compressive stress and impact 



o 



C/3 

O I velocity, sliding occurs in a crack- like mode. For higher values of the initial com- 

c/3 ' pressive stress and/or impact velocity, sliding takes place in a pulse-like mode. A 

^~». variety of sliding modes are obtained in the calculations depending on the impact 

o ' velocity, the initial compressive stress and the values of interface variables. One 

pulse-like mode involves well-separated pulses with the pulse amplitude increasing 
with propagation distance. Another pulse-like mode involves a pulse train of essen- 
tially constant amplitude. The propagation speed of the leading pulse (or of the 
cd ' tip of the crack-like sliding region) is near the longitudinal wave speed and never 

less than \/2 times the shear wave speed. Supersonic trailing pulses are seen both 
experimentally and computationally. The trends in the calculations are compared 
with those seen in the experiments. 
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1 Introduction 

Frictional sliding along an interface between two rapidly deforming solids is a basic prob- 
lem of mechanics that arises in a variety of contexts including moving machinery surface 
interaction (both macro and micro machines), material processing (e.g. cutting), the 
failure of fiber reinforced composites (e.g. fiber pullout) and earthquake dynamics (fault 
rupture). However, a framework for quantifying the wide range of observed dynamic 
frictional phenomena is only beginning to emerge. The classical Amontons-Coulomb de- 
scription of friction states that the shear stress at an interface is proportional to the 
normal stress, with the coefficient of proportionality being the coefficient of friction. Two 
coefficients of friction are identified; a static coefficient of friction that governs the onset of 
sliding and a dynamic coefficient of friction that characterizes the behavior during sliding. 

At the microscale, an evolving population of deforming and fracturing contacts, pos- 
sible phase transitions and the presence of various lubricants play an important role in 
setting the static and dynamic coefficients of friction as well as in governing the transition 
between them. Rate and state models of friction aim at incorporating the effects of these 
microscale processes through appropriately chosen state variables, e.g. Dieterich (1979), 
Ruina (1983), Rice and Ruina (1983), Linker and Dieterich (1992), Prakash and Clifton 
(1993), and Prakash (1998). 

Rate and state models of friction have come to the fore because they substantially 
influence the predicted mode and stability of sliding. Of particular interest is whether 
sliding occurs in a crack-like mode in which the surfaces behind the leading edge of sliding 
continuously slide or in a pulse-like mode, flrst proposed by Heaton (1990), in which sliding 
occurs over a relatively small propagating region. One signiflcance of the sliding mode is 
that the calculated frictional dissipation in the pulse-like mode is signiflcantly less than 
in the crack-like mode and is consistent with some values of heat generation inferred from 
geophysical fleld measurements (Heaton, 1990). 

The classical Amontons-Coulomb description of friction is inadequate for address- 
ing fundamental issues of sliding along interfaces between elastic solids because with 
Amontons-Coulomb friction sliding along such an interface is unstable to periodic per- 
turbations for a wide range of friction coefficients and material properties, with a growth 
rate proportional to the wave number (Renardy, 1992; Adams, 1995). When generalized 
Rayleigh waves exist, Ranjith and Rice (2001) found that there are unstable modes for 
all values of the friction coefficient. Mathematically, instability of periodic perturbations 
renders the response of a material interface with Amontons-Coulomb friction ill-posed. 
Physically, it implies that during sliding energy is transferred to shorter wave lengths, 
leading to pulse sharpening and splitting. Numerically, the splitting of individual pulses 
creates an inherent grid-size dependence (Andrews and Ben Zion, 1997; Ben Zion and 
Andrews, 1998). Ranjith and Rice (2001) showed that an experimentally based rate 
and state dependent friction law (Prakash-Clifton, 1993; Prakash, 1998), in which the 
shear strength in response to an abrupt change in normal stress evolves continuously with 



time, regularizes the problem. Convergence through grid size reduction is then achieved 
(Cochard and Rice, 2000). 

Although rate- and state-dependent frictional laws provide regularization, slip pulses 
may still grow with time. Perrin et al. (1995) and Zheng and Rice (1998) found from anti- 
plane strain (mode-Ill) calculations conditions on the frictional law and the loading con- 
ditions needed for well-posedness. They also found that under certain conditions, which 
include the rate of steady state velocity weakening, sliding can be partially crack-like and 
partially pulse-like. Ben-Zion and Huang (2002) numerically analyzed a configuration of 
two identical solids separated by a variable-width fault zone layer and assumed constant 
external compressive and shear fields. They showed that the friction law used by Ranjith 
and Rice (2001) and Cochard and Rice (2000) regularizes the Adams instability for iden- 
tical as well as for dissimilar materials. The same conclusion was reached analytically by 
Rice et al. (2001). The use of rate and state dependent friction not only regularizes the 
ill-posedness of the sliding problem but can also eliminate supersonic propagation (Rice 
et al., 2001). However, the physical reason for excluding supersonic propagation is not 
clear. 

The issue of limiting propagation speeds arises in the dynamic fracture mechanics 
of growing shear cracks (Freund, 1998; Broberg, 1999; Rosakis, 2002) which has many 
similarities with the frictional sliding process. Early analytical and numerical studies of 
intersonic shear crack growth using continuum slip weakening (Burridge, 1973; Andrews, 
1976; Broberg, 1995, 1996) and velocity weakening cohesive zone models (Samudrala et 
al., 2002) concentrated on elucidating the mechanism of transition from sub-Rayleigh to 
super-shear wave propagation speeds, see also Madariaga and Olsen (2000) and Dunham 
et al. (2003), and on the identification of stable intersonic propagation speed regimes. 
Of particular relevance to the present study is the persistent occurrence of intersonic 
shear rupture along interfaces between identical solids which was first observed experi- 
mentally by Rosakis et al. (1999) and by Coker and Rosakis (2001). This phenomenon 
was subsequently modeled numerically by Needleman (1999) and by Hao et al. (2004). 
Atomistic models of intersonic shear rupture (Abraham and Gao, 2000; Abraham, 2001; 
Gao et al., 2001), field observations of intersonic rupture events during recent large crustal 
earthquakes (Archuleta, 1984; Olsen et al., 1997; Hernandez et al., 1999; Bouchon et al., 
2001; Lin et al., 2002; Bouchon and Valle, 2003), as well as recent laboratory models 
of earthquake rupture (Xia et al., 2004) have demonstrated the remarkable length scale 
persistence (over eleven orders of magnitude) of intersonic rupture phenomena and of the 
main features observed in laboratory experiments and in continuum theories. 

Here, frictional sliding between identical rectangular plates subject to an initial com- 
pressive stress and impact shear loading is studied both experimentally and numerically. 
In the experiments, a uniform compressive stress is applied to two Homalite plates and 
the impact loading is imposed using a gas gun and a steel projectile. Dynamic photoelas- 
ticity is used to record fringe patterns on a micro-second time-scale. The fringe patterns 
give indirect evidence of the shding mode and the propagation speed of the sliding tip 



is measured. Plane stress calculations modeling the experimental configuration are car- 
ried out with friction characterized by a rate- and state-dependent relation, including the 
Prakash-Clifton (Prakash and Clifton, 1993; Prakash, 1998) normal stress dependence. A 
parameter study is undertaken to predict possible modes of sliding in this configuration 
together with the associated propagation speeds. The dependence on loading conditions is 
explored and the computational results are compared with the experimental observations. 

2 Formulation 

2.1 Boundary value problem 

We consider two solids undergoing relative motion across an interface Sint- A framework 
is used where two constitutive relations are specified; one for the plate material and one 
to characterize frictional sliding in terms of the traction and velocity jump across the 
interface. Plane stress conditions are assumed, geometry changes are neglected^ and the 
principle of virtual work is written as 

/ o- : 5edA -I T ■ 5^udS = [ T ■ 6udS - / p— ^ • 5udA (1) 

Ja Js,„t JSext Ja ot^ 

where t is time, cr is the stress tensor, e is the strain tensor, u is the displacement vector, 
T is the traction vector, Au is the displacement jump across the interface, cr : e denotes 
(Tijeji, and A, Sext and Sint are the area, external boundary, and interface line, respectively, 
in the reference configuration. 

Computations are carried out for the specimen geometry shown in Fig. Q] with i = 75 
mm and w = 137 mm. A Cartesian coordinate system with the Xi — X2— plane being 
the plane of deformation is used and the origin is taken as shown in Fig. ^ The plates 
are subject to a uniform compressive stress of magnitude Sq. For t < 0, ^22 = — Sq and 
Tn = So along the interface. Displacements are measured from the uniformly compressed 
state at t = 0~. 

At t = 0^, a normal velocity is prescribed along the portion of the edge Xi = for 
which —h<X2< 0, with 6 = 25 mm, and the shear traction is taken to vanish there. 
Hence, 

ui = - I V{C)di , T2 = on xi = and -6 < X2 < (2) 
h 

^The calculations are carried out using a finite element code developed for finite deformations in 

which the constitutive relation is expressed in terms of second Piola-KirchhofF stress and Lagrangian 

strain components; otherwise finite deformation effects are neglected in the calculations here. 
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\/(f\ } Mmpj tor t-r S t- S t-p /q\ 

^^ ~ ^ Mmp [1 - (t - tp)/t.], for tp < t < (tp + tj ^'') 

0, for t > (tp + ts) 

Here, t^- is the rise time, tp is the pulse time and tg is the step down time. In the 
calculations tr and tg are fixed at 10 /is and (tp — t,.) is 50 /is. On the remaining external 
surfaces of the specimen T = 0. 

The plate material is characterized as an isotropic elastic solid with Young's modulus E 
and Poisson's ratio z/, with properties representative of Homalite: E = 5.2 GPa, u = 0.34 
and density p = 1230 kg/m^. These give elastic wave speeds of 

Q = 2187m/s Cs = 1255m/s CR = 1155m/s (4) 

where q is the longitudinal wave speed, Cs is the shear wave speed and cr is the Rayleigh 
wave speed. 

A rate- and state-dependent friction law is taken to characterize the response along 
the interface X2 = 0. With n and s, respectively, denoting in-plane directions normal and 
tangential to the interface, Aiin and Aug are the difference in velocity between material 
points that are on opposite sides of the interface in the initial configuration. Since the 
change in the positions of material points along the interface is not accounted for, the 
formulation is restricted to small amounts of sliding. When Am„ > — So/C„, the interface 
is in tension and 

T„ = Ts = Q (5) 

along Sint- Otherwise, 

Tn = -CrAUn (6) 

t = Cs[Aus-sgn{Ts)Au,up]. (7) 

Here, as in Povirk and Needleman (1993), the interface is presumed to have an elastic 
stiffness and AusUp is the magnitude of the sliding velocity specified by the rate- and 
state-dependent friction law (AwsHp > 0). In the elastic- friction relation eq. ((7j) AwsHp is 
an internal variable determined from the frictional constitutive description, while Am^ is 
the jump in particle velocity across the interface. Unless otherwise specified, the values 
Cn = 0.3 MPa/m and Cg = 0.1 MPa/m are used in the calculations here. 

The finite element discretization is based on linear displacement triangular elements 
that are arranged in a 'crossed-triangle' quadrilateral pattern. When the finite element 
discretization of the displacement field is substituted into eq. (P) and the integrations are 
carried out, equations of the form 

92 U 

M-^ = R (8) 



are obtained where U is the vector of nodal displacements, M is the mass matrix and R is 
the nodal force vector consisting of contributions from the area elements and the interface. 
Four point Gaussian integration is used along the interface. These finite element equations 
are integrated numerically by an explicit integration procedure, the Newmark /3-method 
(Belytschko et al., 1976). A lumped mass matrix is used instead of the consistent mass 
matrix, since this has been found preferable for explicit time integration procedures, from 
the point of view of accuracy as well as computational efficiency (Krieg and Key, 1973). 
The mesh used in the calculations consists of a fine region with uniform rectangular 
elements near the interface at impact edge, with the mesh spacing then gradually increas- 
ing in size. The mesh has 18, 320 quadrilateral elements and 56, 580 degrees of freedom, 
with the uniform mesh region consisting of 200 x 20 rectangles, each of which is 0.30 
mm X 0.20 mm, so that the uniform mesh extends 60 mm from the impact edge. The 
configuration analyzed and the finite element mesh used are shown in Fig. ^ 

2.2 Rate- and state-dependent friction law 

The classical Amontons- Coulomb frictional relation has the form 

|T,|=/iT„ (9) 

where T„ is the component of T in the direction of the normal to the interface, Tg is 
the component of T tangent to the sliding direction and fi is the (constant) coefficient of 
friction. The relation holds when the normal traction is compressive, i. e. with the 
sign convention here when T^ > 0. 

The problem of two elastic half-spaces sliding with a constant coefficient of friction, 
/i, is ill-posed for a significant range of values of /^, Adams (1995). This manifests itself in 
numerical solutions through a lack of convergence with increasing mesh refinement. Rate 
and state dependent friction laws were introduced, mainly in the geophysics literature, to 
account for experimental observations that could not be rationalized in the context of a 
Coulomb friction description. At each point on the bounding surface, the coefficient of 
friction fi in eq. is taken to depend on the frictional sliding velocity at that point, 
Ausiip, and a set of state variables, di so that 

\Ts\ = fi{6i,Ausiip)Tn 

The state variables are regarded as phenomenological parameters that account for the 
change of contact quality between the surfaces over time. Rate and state dependent fric- 
tion models can account for the following fundamental observations on friction (Dieterich, 
1979; Rice and Ruina, 1983; Ruina, 1983): 

(i) there is an instantaneous increase in the coefficient of friction in response to a step 
increase AiigUp (the direct effect), 
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(ii) there is a subsequent change to a steady-state value of the coefficient of friction with 
this steady-state value being a decreasing function of A-UsUp, and 

(iii) the approach to this steady-state value occurs over a characteristic distance that is 
independent of Augiip- 

Another key observation is that transmitted shear stresses do not instantaneously 
follow a step drop in the normal stresses, Prakash and Clifton (1993), Prakash(1998). 
Instead, immediately following a step drop in |Tn|, as the slip accumulates the shear 
stress gradually approaches a new steady-state level characteristic of the new normal 
pressure and the current sliding velocity, Prakash and Clifton (1993). The relation, eq. 
©, is replaced by, Prakash and Clifton (1993), Prakash(1998), 



|T,|=/i(^o,AM,up)(0i + 02) 
with the evolution of the state variables 9i and 62 given by 

Oi = -^[^i-CT„]Am,hp, 

02 = -^[e2-DTn]Au 



slip- 



(10) 



:iii 



(12) 



where Li and L2 are characteristic lengths. In the steady state, when ^i = and 62 = 0, 
eqs. (Uni) to (O imply T, oc T„. 

In Dietrich (1979), Rice and Ruina (1983) and Ruina (1983) fi was taken to have a 
logarithmic dependence on the sliding velocity. In such a relation, /i is not bounded at 
Awsiip = which leads to difficulties in numerical computations. Rice and Ben-Zion (1996) 
and Ben-Zion and Rice (1997) appealed to an activated rate process interpretation which 
resulted in a relation bounded at AwsUp = 0. Here, we adopt a purely phenomenological 
expression for fi proposed by Povirk and Needleman (1993) that has the form 



/i(6'o,Au 



slip J 



jwf^ + l 



1/m 



Vn 
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where the evolution of 6*0 is given by, 

60 = 3(1 



^oA-itsiip 



(15) 



and fig, fJ'd, Lq, Vq, V\ and B are prescribed constants. 

The steady state value of 6*0, i.e. the value at which ^q = is Lo/A-UsHp. Substituting 
this value into eqs. (|T3|) and (|TH) gives the steady state value for the friction coefficient, 
[iss-, at sliding velocity A-UsUp as 



[i'ss = /id + {^s - ^d) exp 






(16) 



The frictional relation eq. (|7j) is only invoked if there is contact (T„ > 0). Even when 
there is contact, AugUp may still vanish. Define 

then from eqs. ([111)) and (fT^ . and requiring AwsHp to be non- negative, AugUp is given by 

_fK,(/3"^-l) for/3>l 
A^siip - I Q for /? < 1 ^^^) 

Note that /? > 1 implies \Ts\ — g{6o){6i + O2) > 0. When the normal traction is constant 
9i + 02 = Tn and /? simplifies to \Ts\/Tng{9o) as in Povirk and Needleman (1993). 

We refer to circumstances where AwsUp = as sticking and circumstances where 
Ausiip > as frictional sliding. 

The procedure used to carry out time integration of this frictional constitutive relation 
is described in the Appendix. 

2.3 Friction law parameters 

Samudrala et al. (2002) and Rosakis (2002) fit sliding experiments of Homalite on Homa- 
lite (Rosakis et al., 1999; 2000) with an expression of the form 



(J'ss — /^o 



Awsiip G 
1 + a — 



2ro 



(19) 



where 2G = E/{1 + u) and tq is the shear strength of Homalite. From the observed 
inclination of tensile micro-cracks emanating from the interface, the parameter values 
a = —0.4, G/tq = 136, and /iq = 0.6 were obtained. 

Parameter values in eq. ()16|1 were chosen to match the response in eq. (I19j) for A-UsHp < 
40 m/s. Figure Et shows fiss as a function of A-UsHp from eq. (|T!l|) and from eq. (fTHI) with 
fis = 0.6, ^d = 0.5, p = 1.2 and Vi = 26 m/s. The parameters Vq and m in eq. (jT^ 
were chosen to be consistent with the observed response for AwgUp ranging from 10 m/s 
to 40 m/s (Fig. ^). Figure Eb shows the friction coefficient, Tg/Tn, as a function of of 
the accumulated frictional sliding AusUp = J AusUpdt. There is an instantaneous increase 



in Tg/Tn in the same direction as the change in AwsUp (the direct effect) followed by a 
gradual decrease to a steady-state value. The apparent coefficient of friction TgjT^ can 
attain values outside the range defined by the static and dynamic coefficients of friction 
due to the direct effect in the rate- and state-dependent friction law. 

The parameter values for the normal stress dependent response, C, D, Li, are chosen 
close to values obtained in the experiments of Prakash and Clifton (1993). The initial 
values for the internal variables were set to ^o(O) = L/lOVi and 6*1(0) = CSq and ^2(0) = 
DT,Q. Figure El shows the response using eqs. (fTUI) to (fT^ under a step jump in compressive 
normal traction at a constant sliding velocity. The plots show the imposed normal traction 
(dashed line) and the shear traction (solid line) as functions of the accumulated frictional 
sliding Awsiip. There is a gradual change of the state variables and the shear resistance 
after a step decrease in the normal compressive stress. The apparent coefficient of friction, 
Tg/Tn, however, has a sudden jump. 

The choice of the initial values for the state variables is arbitrary but can be of signif- 
icance when comparing to other friction models such as slip-weakening models (Bizzarri 
et al. 2001). The parameters characterizing the frictional constitutive relation used in the 
calculations are those in Table Q] unless specifically stated otherwise. 

3 Experimental Methods and Results 

The experimental procedures are similar to those used to study shear crack propagation 
(Rosakis et al., 1999). Two Homalite plates, each with £ = 139.7 mm w = 76.2 mm (see 
Fig. ^) and 9.525 mm thick are held together by a uniform compressive stress. Homalite 
is a brittle polyester resin that exhibits stress induced birefringence and is mildly rate 
sensitive. At a strain rate of 10'^ s~^ the elastic wave speeds are given by eq. (^. The 
impact loading is imposed via a cylindrical steel projectile of diameter 25 mm and length 
51 mm fired using a gas gun with impact speeds ranging from 10 m/s to 60 m/s. A steel 
buffer plate, 25.4 mm x 76.2 mm and 9.525 mm thick, is bonded to the impacted plate 
to prevent shattering at the impact side (xi = in Fig. ^) and to induce a more or 
less planar loading wave. A uniform compressive stress is applied by a press which was 
calibrated using a load cell. The loading wave, as measured from a strain gage glued 
to the specimen, is of a trapezoidal form with a rise time of 10-20 /xs followed by an 
essentially steady velocity for 40 /xs. Dynamic photoelasticity is used to extract stress 
field information around the interface. The photoelastic fringe patterns were recorded in 
real time using a high-speed Cordin CCD camera capable of capturing 16 images at a rate 
of 100 million frames per second. A collimated laser beam with a diameter of 130 mm is 
used to illuminate the specimen. Two pairs of circular polarizer plates, placed on either 
side of the specimen as described by Rosakis et al. (1999), produce isochromatic fringes. 
The photoelastic optical setup is arranged for light field. 

The isochromatic fringes are related to contours of ai — (T2, with o"i being the maximum 



in-plane principal stress and cr2 the minimum in-plane principal stress, through the stress 

optical relation 

NF 
a,^a, = — ^, (20) 

h 

where F^ is the stress optical coefficient of Homalite, h is the specimen thickness, and N 

is the isochromatic fringe order. 

3.1 Experimental results 

The experimental results are presented in detail in Lykotrafitis et al. (2004). Here, a few 
experimental results that illustrate characteristic features are presented for comparison 
with the numerical calculations. 

Results with a compressive stress of 9.4 MPa and impact velocities of 32.7 m/s and 
42.2 m/s are presented. Figure |3] shows the isochromatic fringe patterns obtained with 
an impact velocity of 32.7 m/s. Experiments conducted at lower impact velocities exhibit 
similar characteristics. In Fig. EJ the loading wave front, which travels at q, arrives from 
the left. Behind the loading wave front a shear Mach cone, formed by a sharp change in 
fringe density is observed. The inset to Fig. |3] focuses on the region near the propagating 
tip with a thick line drawn over the Mach line. The sliding tip follows the loading wave 
and the tip can be located by tracing the Mach line to the interface. In this case the 
sliding tip follows shortly behind the dilational loading wave. The propagation speed of 
the sliding tip is ~ 1810 m/s as measured by two methods. One measurement involves 
following the position of the sliding tip in various frames and knowing the frame timing, 
the propagation speed is obtained. In the second method the Mach angle is used to obtain 
the tip speed. The two methods give consistent propagation speeds. 

In Fig. m there is a concentration of isochromatic fringes some distance behind the 
sliding tip that propagates at approximately cr and that is a consequence of loading 
wave reflections. Other characteristic features in Fig. |3] are: (i) there is a cusp in the 
stress contours at the interfaces, indicating that the propagation speed is slightly faster 
along the interface than in the bulk; (ii) the fringe density is higher in the plate where 
the impact loading is applied, showing that energy is not transferred easily through the 
interface; and (iii) the fringe discontinuity at the interface shows that sliding is occurring 
in a crack-like mode. 

Isochromatic fringe patterns for an experiment with an impact velocity of 42 m/s are 
shown in Fig. El at 40 /xs, 48 /xs and 60 /xs after impact. The thick lines in the insets are 
drawn to show the Mach line locations and orientations. Both methods of calculating 
the propagation speed of the sliding tip give a propagation speed of about 1950 m/s. 
Although the general characteristics seen in Fig. |3] are preserved, new features enrich the 
picture. Behind the impact wave a shear Mach line emanating from the sliding tip is 
observed. In addition, a second Mach line emanates behind the sliding tip that is not 
parallel to the ffist one (Fig. Eli). This Mach line is at a shallower slope corresponding to 
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a supersonic propagation speed of ~ 2600 m/s. In Fig. |S|d the tip of the second Mach hne 
approaches the tip of the first Mach hne at the interface. These two points merge as the 
second point catches up with the first point (Fig. Efc) and only one Mach hne continues 
to be observed in subsequent frames (not shown). 

Behind the second Mach hne there is a fringe concentration, associated with the load- 
ing, that travels at the Rayleigh wave speed. Some fringes pass continuously from the 
upper plate to the lower plate which is consistent with a contact region forming. In exper- 
iments with the magnitude of the initial compressive stress increased, there is an increase 
in the number of discontinuous regions of fringe concentration that may indicate sepa- 
rate regions of sliding and regions of sticking (i.e. multiple pulses). In each experiment, 
the speed of the leading sliding tip is found to be essentially constant. In addition, the 
propagation speed of the leading sliding tip was found to increase with increasing impact 
speed (or decreasing compressive stress) and in all cases to be between a/2cs and q. 

4 Numerical Results 

The focus here is on illustrating the range of behaviors that are obtained by varying 
the magnitude of the initial compressive stress Sq and the impact velocity Mmp- Results 
are presented for combinations of Sq and l^imp that give rise to the following five sliding 
modes: 

Case I: A crack-like mode - Sq = 6 MPa, I^mp = 2 m/s. 

Case II: A pulse-like mode - Sq = 30 MPa, V^mp = 2 m/s. 

Case III A train of pulses - Sq = 10 MPa, V^mp = 20 m/s. 

Case IV: Multiple pulses coalescing to form a crack - Sq = 0.9 MPa, V^mp = 10 m/s. 

Case V: A pulse-like mode followed by a crack-like mode - Sq = 40 MPa, V^mp = 2 m/s. 

In all cases but Case V, the friction parameters are as specified in Table Q For Case 
V, the exponent p in eq. ()14|) is taken to be 0.5 and V^i = 1 m/s. Also, for Case IV the 
values of the interface elastic constants are C„ = 0.03 MPa/m and Cg = 0.01 MPa/m, 
one tenth their value in all the other cases. In the results presented here, the frictional 
sliding region remains within the fine uniform part of the mesh. 

Results are presented for: (i) propagation speeds; (ii) the spatial distribution and time 
evolution of various interface quantities; and (iii) isochromatic fringe patterns (contours 
of CTi — cr2 with (Ti being the maximum in-plane principal stress and (72 the minimum 
in-plane principal stress). 

Curves of propagation speed versus time for all five cases analyzed are shown in Fig. El 
For reference, the longitudinal, q, shear, c^, and \/2cs (= 1775 m/s) wave speeds are 
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marked. Also shown in the figure are data points corresponding to measured propagation 
speeds for the experiment with Eg = 9.4 MPa and an impact velocity of 42 m/s. The first 
step in computing the propagation speeds is to record the location of the point furthest 
from the impact edge where AugUp > 0. This gives the position of the leading sliding 
tip at various times and the propagation speed, Kip, is then calculated from this data 
by a progressive least squares fit using five points. At least some of the oscillations in 
propagation speed in Fig. El may be due to this numerical differentiation procedure. 

Sliding generally initiates somewhat away from the impact edge and then occurs over 
several small sliding regions that link up. Thus, the process of sliding initiation is not 
one of propagation of a sliding region but involves jumps in the location of the leading 
sliding edge. As a consequence, the 'speeds' calculated during the early stages of sliding 
are unphysically high. Subsequently, the propagation speeds decrease and for Cases II 
and V, the mean V^ip is about q whereas for Case IV, after the initial stages of sliding, 
VJ;ip reduces to ~ 2000 m/s, in good agreement with the experimental value of 1950 m/s. 
There are multiple pulses in Case IV and the tip of the second pulse travels at about 2700 
m/s, as compared with ^ 2600 m/s in the experiments. It is worth emphasizing that both 
the experiments and the calculations indicate that the second pulse travels at a supersonic 
speed. However, for the calculation of Case IV, Eq = 0.9 MPa and V^mp = 10 m/s as 
compared with 9.4 MPa and 42 m/s in the experiments. Hence, there is consistency of 
propagation speeds for this mode of interface sliding although the loading conditions to 
achieve this are quite different in the calculations and the experiments. One possibility 
for this discrepency is that the effect of the impedence mismatch between the steel plate 
bonded to the specimen and Homalite is not accounted for in the computations. 

In the following, distributions are shown both for the evolution of Aiig, which is the 
particle velocity jump, and for AwgUp, which is the sliding rate obtained from the frictional 
constitutive relation. From eq. ^, the difference between these is the 'elastic' sliding rate 
Ts/Cs- There are cases where Aiis ~ AusUp with the elastic sliding rate being negligible. 
However, in other cases the difference between Aiis and Augiip is substantial and the 
elasticity of the interface then plays a significant role. 

The principal stress contours that evolve are a consequence of both the imposed loading 
and sliding along the interface. To provide a perspective on the contribution of interfacial 
sliding, the contours of ai — (72 that emerge from impact loading with no sliding are 
considered first. 

4.1 Symmetric loading with respect to the interface 

In order to illustrate the isochromatic fringe patterns that result without sliding, the 
symmetric loading configuration shown in the inset of Fig. [7^ is analyzed. Figure [7^ 
shows the simulated isochromatic fringe pattern at t = 16 fis for a calculation carried 
out with Eq = 10 MPa and \4mp = 5 m/s. At this time, the loading wave front has 
propagated 34 mm into the specimen. The isochromatic fringe pattern for a symmetric 
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loading experiment with Sq = MPa and Vimp = 58 ni/s is shown in Fig. [7|d. In both 
the experiments and simulations the continuity of fringes across the interface behind the 
loading wave front indicates that no sliding occurs along the interface under symmetric 
loading. 

In both the simulations and the experiments, at the loading wave front there is a 
cusp in the stress contours at the interface (marked by the arrow A), indicating that the 
propagation speed is faster along the interface than in the bulk. The cusp in the fringes at 
the interface in Fig. [7|d is also seen when frictional sliding occurs (Figs. 0] and Ej). Another 
characteristic feature is marked by box B in Fig. [3 The curvature of the fringe lines 
associated with this feature is greater in the computations in Fig. EK, but this curvature 
decreases with time. A similar feature, with a curvature closer to what is seen in the 
computations, is in Fig. EId. Our simulations indicate that the cusp in the fringe line and 
the feature in box B are a consequence of the impact loading conditions and thus not 
directly associated with frictional sliding. 

A relatively broad head wave that emanates from the interface following the loading 
wave front is shown by the arrow C The width of the head wave in the simulations was 
found to depend on C„ and C^ in eq. ((Tj) and thus can be used to set the values of the 
interface elastic stiffnesses. 

4.2 Case I: Crack-like mode 

Figure IHfe' shows the variation of Awsiip, and the shear traction, T^, along the interface at 
t = 22 yus (dashed lines) and at t = 32 fis (solid lines) for Case I (Eq = 6 MPa, V^mp = 2 
m/s). Frictional sliding initiates at the impact edge (xi = 0) 16 /xs after impact. The 
distributions of A-UsUp and T, along the interface are similar at 22 /xs and 32 /is, indicating 
at least quasi-steady behavior. 

Ahead of the frictional sliding tip, AusUp = 0. At the tip, AwsUp jumps to its peak 
value of about 7 m/s within two mesh spacings (eight interface integration points) and 
then gradually reduces to ~ 5 m/s. As seen in Fig. |Hk, the peak value of A-Usiip gradually 
increases with time. This behavior is crack-like in that frictional sliding persists at a point 
on the interface after the sliding front has past that point. 

The shear traction steadily increases from the loading wave front, which is at xi = 71 
mm at t = 32 fis, and reaches a maximum value of about 3.9 MPa at the frictional sliding 
tip after which it decreases to a steady value of 3.6 MPa. The normal traction along the 
interface gradually increases from its initial value of 6 MPa to 6.39 MPa at the frictional 
sliding tip and subsequently decreases to 6.14 MPa along the sliding region. The apparent 
coefficient of friction, /Xapp = Tg/Tn, is ~ 0.609 at the frictional sliding tip and reduces to 
0.585 along the sliding region with small variations. Thus, in this case /iapp remains very 
close to yUg. The value of T„ itself is nearly constant along the interface, differing by less 
than 0.6 MPa from Eg = 6 MPa, except near the impact edge (0 < Xi < 10 mm) where 
the magnitude of T„ decreases. 
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The actual velocity jump across the interface is Aug, which consists of an elastic part 
in addition to A-UgUp- Figure IHb shows the distribution of Aug as well as the traction 
increment dTg = Tgdt distribution. Because of the elasticity in the interface constitutive 
relation, eq. ((Zj), Aiig increases gradually to it peak value, rather than exhibiting a sharp 
jump. However, behind the frictional sliding tip, Aus and AitgUp are nearly equal. Consis- 
tent with this, since the elastic contribution is Tg/Cg, Tg ^ behind the frictional sliding 
tip. 

Isochromatic fringe patterns (contours of cti — (J2) at t = 32 us are shown in Fig. El 
Behind the loading wave front (which is at 71 mm) a concentration of fringes at 44 
mm indicates the location of the frictional sliding tip. A shear Mach cone, formed by a 
sharp change in fringe density, emanates from the sliding tip at an angle of 34°. Using 
Aa/At = sin6/cg, where 6 is the Mach angle, gives a tip speed of 2244 m/s which is 
close to, but less than, the directly calculated tip speed in Fig. El Another Mach cone, 
formed by a more gradual change in fringe density at an angle of 67° and, presuming that 
Q is the relevant speed of sound for this Mach cone, a propagation speed of 2360 m/s is 
obtained from Aa/At = sin6/ci. Behind the frictional sliding tip there is no significant 
stress concentration, indicative of sliding occurring along the interface. 

4.3 Case II: Pulse-like mode 

In Fig. ITUk . the variation of the frictional sliding rate, Aiisiip, and the shear traction, Tg, 
along the interface at t = 32 fj,s (dashed lines) and at t = 34 /is (solid lines with circles) 
are shown for Case II (So = 30 MPa, Vimp = 2 m/s). Frictional sliding initiates at t = 28 
/is, which is later than for Case I because the larger compressive stress requires a larger 
shear stress for initiation. In Case II, frictional sliding occurs in a pulse-like mode which 
involves a relatively narrow zone (the pulse) where AiisUp > 0, with A-itsUp = ahead of 
and behind the pulse. The circles in Fig. ITUk show the integration point locations (since 
Ausiip is a constitutive quantity it is stored at interface element integration points) to 
indicate the extent to which the pulse is resolved; the sharp pulse consists of 9-10 points. 

A pulse initiates at the impact edge and as it travels across the interface the peak 
value of Awsiip increases. The pulse profile consists of AusUp gradually increasing from 
zero to 40 m/s in 4.5 mm after which it sharply rises to 700 m/s over a distance of 0.5 
mm and then drops abruptly to zero over a distance of 0.4 mm. Additional pulses initiate 
and evolve in a similar manner; there are two pulses at t = 32 /is and three at t = 34 /is. 

The shear traction distribution along the interface is also shown in Fig. ITUk (note that 
the range for Tg is larger than in Fig. |S^). At t = 32 /is, the shear traction along the 
interface gradually increases from zero at the loading wave front (which is at 71 mm) to 18 
MPa at the tip of the leading pulse. The value of Tg remains approximately constant over 
the pulse, abruptly decreases to 14 MPa at the trailing edge of the pulse and then very 
gradually increases until the second slip-pulse is reached. The variation of the apparent 
coefficient of friction, yU^pp = Tg/T^, is quite different from that for the crack-like mode. 
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Here, at t = 32 /is, /iapp is 0.60 at the front of the frictional shding tip, very shghtly 
increases to 0.61 over the pulse and then, with the drop in Tg, falls to 0.44. In this 
case, neither fig nor fid serve as bounds on the apparent coefficient of friction. Behind 
the leading pulse the apparent coefficient of friction steadily increases to 0.60 until the 
trailing pulse is reached. The magnitude of T„ is nearly constant at Sq = 30 MPa over 
most of the interface but decreases near the impact edge (0 < Xi < 10 mm). 

In contrast to the crack-like mode, for the pulse-like mode there is a qualitative differ- 
ence between the distributions of AiisUp and Aug along the interface. As seen in Fig. ITUb . 
which shows distributions of Aiis and dTg = Tgdt at t = 34 /is, no pulse occurs in the 
distribution of Aug- However, there is an abrupt change in Aiig at the location of each 
pulse. On the other hand, there are pulses in the traction increment dTg and hence in the 
elastic contribution to Aug = Tg/Cg. The abrupt change in Aiig corresponds to a pulse 
in dAug/dx, i.e. to a 'weak' pulse as opposed to the 'strong' pulse in A-UsUp- 

Isochromatic fringe patterns (contours of cxi — (72) at t = 32 fis are shown in Fig. [TT] 
where only the stress contours behind the loading wave front are visible in the region 
shown. There is a concentration of fringes at approximately Xi = 22 mm, which corre- 
sponds to the location of the leading slip pulse. The change in the fringe spacing and 
concentration along a shear Mach line of 26° indicates a propagation speed of about 
2800 m/s (sin 6'/cs). The dense set of fringe near the impact site arises from a slip pulse 
initiating there. 

4.4 Case III: A train of pulses 

Figure IT2k shows distributions of frictional sliding rate, A-UsUp, and shear traction Tg for 
a train of pulses (Sq = 10 MPa, \4mp = 20 m/s). The distributions are shown at two 
times, at t = 32 fis and at t = 33 /is. As in Fig. ITUk . the circles indicate the positions 
of the integration points at which Ausiip and Tg are evaluated. A key difference between 
the pulses in Case HI and those in Case II is that in Case HI the peak value of Awsiip, 
which is ~ 120 m/s, does not increase with propagation distance or time. The variations 
in the peak value of AiisUp seen in Fig. IT^ is mainly due to some pulses attaining their 
peak value between the plotted integration points. Typically, a pulse extends over four 
elements (sixteen integration points). The slip pulses at 33 /is (dashed lines) lie on top 
of the pulses at 32 /is with a one pulse offset, illustrating the self-similar propagation of 
these pulses. 

The shear traction along the interface has a mean value of about 7 MPa with oscil- 
lations of 1.5 MPa. The apparent coefficient of friction, /iapp = Tg/Tn oscillates between 
0.59 and 0.48 for the slip-pulses behind the leading pulse, except for the leading pulse 
where the maximum is 0.63. The value of T„ is fairly uniform over the interval shown, 
having a maximum of about 14 MPa at xi ~ 45 mm and decreasing smoothly to around 
11.5 MPa at Xi = 30 mm and Xi = 60 mm, so that the variation in /iapp arises principally 
from the oscillations in T,. Due to the direct effect in the friction relation, there is a 
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lag between the onset and cessation of frictional sliding and the maximum and minimum 
values of /iapp- The maximum value of 0.59 for /iapp is reached 0.43 mm behind the sliding 
tip for the pulses following the leading pulse and the minimum value of 0.48 occurs 0.25 
mm behind the location where AusUp peaks (for all pulses, including the leading pulse). 

As also seen in Fig. [T2k, except for the leading pulse, the shape of the pulses is the 
same, where the width of the frictional sliding region (A-UgHp > 0) is 1.33 mm and the 
width of the region with AwsUp = is 0.92 mm. The spacing between the pulses is 2.21 
mm. The leading pulse has more structure to its front edge, and this structure is preserved 
as the pulse propagates. 

The velocity jump across the interface, Aiis, and the traction increment dTg = Tgdt 
are shown in Fig. IT^ . The traction increment distribution reflects the pulse structure of 
Awsiip; while the oscillations in A-itg are smoother. 

Figure IT^ shows the isochromatic fringe patterns (contours of 0"i — (12) at t = 32 /is. 
At the front of the first pulse at 55 mm, a shear Mach wave occurs at 37°, corresponding 
to a propagation speed of ~ q. Behind the slip front vortex like stress contours are seen 
that correspond to the back and front of each pulse together with shear Mach waves at 
each transition from slip to stick. Behind the propagating front, there is a vortex-like 
structure. The change in the contours at X2 = 2 mm is due to the change in resolution 
of field quantities associated with the transition from the uniform mesh to the coarser, 
graduated mesh. This calculation was also carried out for a finer mesh in which the 
uniform mesh region extends to X2 = 4 mm, with the element size in the uniform region 
fixed. The results are essentially unchanged, but the finer mesh results in the details of 
the stress distribution behind the propagating front being visible over the larger region 
resolved. 

The accumulated frictional sliding, AwgHp = /q A-UsHprft, the displacement jump across 
the interface, Aug and the shear stress T^ distributions for the crack-like propagation 
mode (Case I) and for the train of pulses (Case III) are compared in Fig. El In both 
cases, the difference between Aug and AwsUp is due to the assumed interface elasticity. In 
Fig. [T3K for Case I, both Aug and AtisUp vary smoothly along the slipped surface. For 
Case III, Fig. IT^ . although the displacement jump itself varies smoothly, steps in AwsUp 
are evident. The elasticity smooths out the steps in Am^ along the interface, but these 
are reflected in the oscillatory distribution of shear traction, T^, along the interface. 

4.5 Case IV: Multiple pulses coalescing to form a crack 

In Case IV, the values of the interface elastic constants are C„ = 0.03 MPa/m and 
Cs = 0.01 MPa/m, one tenth their value in all the other cases. Figure ITSl shows the slip 
velocity A-UsUp along the interface at 22 /xs, 32 fis, and 41 /xs. Sliding initiates at 14 /is 
and at 22 /is (Fig. IT^ ) three pulses are observed. The tip of the leading pulse at 18 mm 
and the distance between the end of the first pulse and the front of the second pulse is 3.4 
mm. As in Cases II and III, the leading pulse includes a gradual rising part in the front 
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before there is a steep rise in A-UsUp whereas the second and third pulses do not exhibit 
such a gradual rise in A-UgUp- At 32 /is in Fig. HSb, the tip of the leading pulse has moved 
to 41 mm with the distance between the first and second pulse decreasing to 0.4 mm. 
During this time the average speed of the tip of the leading pulse is ~ q and the average 
speed of the second pulse is supersonic, ~ 2700 m/s. Thus, the second and third sliding 
pulses which are traveling at supersonic speeds catch up with the lead pulse and coalesce 
to form a single crack-like sliding region as shown in Fig.fT^b (at 41 /is) and then continue 
to travel at q. This crack-like region expands with time as more pulses are generated at 
the impact edge, travel at supersonic speeds and coalesce with the main sliding region. 

Isochromatic fringe patterns (contours of 0"i — (J2) are shown in Fig. ITHk . b, and c, at 
22 yus, 32 /is, and 41 /is, respectively. The contours are shown for X2 > because the 
Mach line are clearer in this plate. At t = 22/is the tip of the sliding pulses is at 18 mm 
as shown by the change in the spacing of the stress contours in Fig. ITBk . In Fig. ITBb . 
the tip has moved to 41 mm. There are three regions where an abrupt change in the 
density of the contours takes place forming three shear Mach lines at three angles to the 
interface as shown in the inset. The intersections of these lines with the interface are at 
the locations of the tips of the three sliding pulses shown in Fig. IT^ . The shear Mach 
angles of 34°, 29.3° and 23.4° correspond to propagation speeds of 2190 m/s, 2560 m/s 
and 3160 m/s for the first, second and third pulses, respectively. In Fig. IT6h . at 41 /is, 
the stress contours show a concentration of fringes only at 60 mm corresponding to the 
sliding tip. At this time the second and third pulses have coalesced with the first pulse to 
form a crack-like sliding region propagating at the longitudinal wave speed. For xi < 25 
mm there are disturbances that are pulses traveling at supersonic speeds that eventually 
coalesce with the crack-like sliding region, thereby expanding it. 

4.6 Case V: A pulse-like mode followed by a crack-like mode 

In Case V, p = 0.5 and Vi = 1.0 m/s in eqs. (fT^ and (fTB|) . Sliding initiates as a slip pulse 
at 31 /is. Figure [121 shows the distribution of AiisUp along the interface at three times; 
t = 32 /is, t = 38 /is and t = 44.6 /zs. Also shown is the distribution of the shear traction, 
Ts, along the interface at t = 44 /is. A growing pulse is obtained with the pulse width 
decreasing from 8.5 mm at t = 32 fis to 2.5 mm at t = 44 /is and the peak value of AitsUp 
increases from 155 m/s to 392 m/s. The propagation speed of the pulse varies slightly 
over the range of times in Fig. El 2600 m/s at t = 32 /is and 2450 m/s at both t = 38 
/is and t = 44.6 /xs. After some oscillation behind the growing pulse, AitsUp is constant at 
about 21 m/s. The shear traction increases to a maximum value of 23.5 MPa at the tip 
of the pulse, then drops to 18.5 MPa in the pulse region, after which it remains almost 
constant at 20 MPa over the rest of the sliding region. The apparent coefficient of friction 
changes from a maximum value of 0.595 at the sliding tip to a minimum value of 0.47 at 
the point of maximum frictional sliding rate until it settles at a steady state value of 0.50 
along the rest of the sliding surface. 
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The isochromatic fringe patterns (contours of o"i — (T2) are shown in Fig. ^Jat three 
times of 32 /xs, 38 /xs and 44.6 /is. The arrow in each figure point to a hne (drawn off the 
interface) that indicates the width of the pulse and its position along the interface. The 
pulse width decreases as the crack catches up to the pulse with time and distance. At 
t = 38/is (Fig. CHb), there are two separate Mach waves corresponding to the slip pulse 
and the initiation of crack-like slip. As the pulse width diminishes in length the tip of the 
crack-like region approaches the leading tip of the pulse. From the shear Mach angles, 
the propagation speed of the back of the pulse in Fig. ITHb is approximately 3500 m/s, 
while the front of the pulse is traveling at about 2600 m/s. 

4.7 Dependence of the sliding mode on the initial compressive 
stress and impact velocity 

The dependence of the frictional sliding mode on the initial compressive stress, Eq, and 
the impact velocity, V^mp, is shown in Fig. UHl All calculations in Fig. IT^ were carried out 
using the interface friction properties listed in Tableland C„ = 0.3 MPa/m, C^ = 0.1 
MPa/m. The dashed line indicates the boundary between the occurrence of crack-like and 
pulse-like modes. It is near the dashed line that mixed crack-like and pulse-like modes 
occur. Quite generally, for low impact velocities and low values of the applied compressive 
stress, the crack-like mode is obtained. On the other hand, for a high magnitude of the 
initial compressive stress and a sufficiently low impact velocity, the growing pulse mode 
occurs. 

Figure EDI illustrates the nature of the transition for a fixed value of Sq and varying 
l^mp, corresponding to the line A in Fig. UHl In Fig. 1201 ^0 is fixed at 10 MPa and \4mp 
varies between 2 m/s to 40 m/s. Values of Aiis, A-UsUp, Ts and /Xapp = Ts/Tn at xi = 16.1 
mm on the interface are plotted as functions of time. With V^mp = 2 m/s, essentially 
steady values of the velocity jump across the interface, Aiig, and the frictional sliding 
rate, Ausnp, are attained as seen in Figs. EUk and OUb . respectively. With the impact 
velocity increased to 10 m/s, the mode has changed to one involving a train of pulses, as 
seen in Fig. 120b. With a further increase in impact velocity, the pulse amplitude increases 
as does the pulse frequency. Correspondingly, the oscillations in the velocity jump across 
the interface increase. However, as seen by comparing the response for ^mp = 20 m/s 
with \^mp = 40 m/s in Fig. I^Hk . the amplitude of the oscillations in Aiig do not continue 
to increase with increasing impact velocity, in contrast to the pulse amplitude in Fig. 120b. 
Rather, there is an increase in the 'elastic' part of the velocity jump and this is refiected 
in the increased magnitude of the shear traction Tg in Fig. I^Ub . The oscillations in the 
apparent coefficient of friction, yUapp = Ts/Tn, increase in magnitude with increasing V^mp! 
yUapp changes from an almost steady value of 0.598 at \4mp = 2 m/s to /iapp oscillating 
between 0.598 and 0.470 at l^mp = 40 m/s. In addition, the frequency of the oscillations 
in /iapp increases with increasing impact velocity. 

The nature of the transition for V^mp fixed at 2 m/s and with varying Sq, the line 
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B in Fig. ^J is illustrated in Fig. |^ As in Fig. 1201 the evolution with time of Aiig, 
Awsiip, Tg and /^app at a fixed point along the interface (xi = 16.1 mm) is shown. With 
increasing initial compressive stress, frictional sliding initiates later (Figs. l2Tk and 121b). 
For both Sq = 6 MPa and Sq = 10 MPa, crack-like behavior is obtained. The values of 
Aits and AugUp at the beginning of frictional sliding exceed their subsequent steady-state 
values, which is a consequence of the direct effect in the rate- and state-dependent friction 
relation. With Sq = 16 MPa, a train of pulses in frictional sliding rate (Fig. I^Tb) occurs, 
with oscillations both in Aug fFig. I^Tk) and in Tg fFig. I?lh). When the compressive load 
is increased to Sq = 30 MPa, steeper and narrower pulses are obtained. Fig. I^Tb . and the 
peak value of AiisUp reaches ~ 820 m/s, much exceeding the scale in Fig. I^Tb . These slip 
pulses increase in amplitude as they propagate and the peak AugUp can be of the order 
of Cs- For a sufficiently high value of Sq, the closely trailing slip pulses are suppressed. 
Then, a trailing slip pulse only occurs after some distance and time has passed so that the 
pulses are well-separated. Such growing pulses have been obtained in studies of sliding 
between dissimilar elastic solids using rate- and state-dependent friction laws (Ben-Zion 
and Huang, 2002; Adda-Bedia and Ben Amar, 2003). It is worth noting that if the scale 
of observation is sufficiently small compared with the pulse spacing, only a single pulse 
is observed which resembles the slip pulses discussed in relation to earthquakes (Heaton, 
1990). 

The nature of the slip pulses that occur when Sq is increased (line B in Fig. ITT?|l differ 
from those obtained when V^imp is increased (line A in Fig. E)). With fixed Sq, increasing 
V^mp increases the magnitude of the slip pulses while the pulse spacing and the pulse 
width decrease. For Eq = 16 MPa, it was found that the pulse width to pulse spacing 
ratio remained constant with varying Vimp and this ratio was equal to one-half. This 
behavior is analogous to the oscillatory stick-slip behavior observed by Baumberger et al. 
(1994). The extent to which the pulse width to pulse spacing ratio depends on Sq remains 
to be determined. In both Figs. OUh andl^Tb. the peak values of the apparent coefficient 
of friction do not increase with increasing Vimp (Fig- llOt) or Sq (Fig. HTh ) even though 
the magnitude of Tg does, thus implying that there is a corresponding increase in T^. 

To illustrate the evolution of sliding in a transition mode, we describe this evolution 
for the calculation with So = 10 MPa, \4mp = 7 m/s (marked by one of the triangles in 
Fig. IT^. Oscillations occur in A-UsUp along the interface, with the highest peak near the 
sliding tip and decreasing in amplitude behind the tip. When sliding has progressed about 
50 mm along the interface, the amplitude of the first three oscillations is large enough 
for Awsiip to vanish at the crests leading to sticking regions (however the maximum value 
of Aiisiip generally remains a small fraction of an elastic wave speed). The amplitude of 
the following oscillations is much smaller and sticking does not occur. Thus, we obtain a 
crack-like sliding region following three slip-pulses of decreasing magnitude. 

Other loading combinations of Eq and Vimp on the transition line in Fig. ^J lead to 
similar oscillations in the slip rate or to unsteady behaviors, such as mixed crack-like 
and pulse-like sliding modes that evolve with both time and distance. This behavior 

19 



is reminiscent of that found in the analysis of spring-shder system obeying a rate- and 
state-dependent friction law (Ruina, 1983; Gu et al., 1984) where it was shown that 
the transition from steady sliding to stick-slip goes through several oscillations, period 
doubling and self-sustained periodic oscillations. Gu and Wong (1991,1994) also found 
this transition to be dependent on the load point driving velocity in addition to the 
stiffness and constitutive parameters. 

5 Discussion 

The experiments and calculations for symmetric loading with no sliding, permit the fea- 
tures of the response associated with sliding to be identified. Also, the good agreement 
between the experiments and the calculations for symmetric loading provides confidence 
that the wave propagation aspects of the experiments are well-represented in the calcu- 
lations. Features that are common between the symmetric loading experiments and the 
corresponding calculations include the cusp in the photoelastic fringes at the interface, 
the formation of a head wave emanating from the loading wave front and the distortions 
in the photoelastic stress field behind the loading wave due to the trapezoidal nature of 
the input. 

A guideline for assessing mesh resolution (Zheng and Rice, 1998) is that the mesh 
spacing (0.3 mm along the interface in the uniform mesh region) should be significantly 
less than a characteristic length, h*, associated with the friction law, 

h* = \2GLo/[7rAu^iip{dTss/dAu^iip)]\ 

where G is the shear modulus and Tgg is the shear stress under steady state sliding. For 
a normal stress of 10 MPa, h* attains its minimum of 0.81 mm for AusUp = ^i = 26 m/s. 
Using Ausiip = V^i = 7 m/s as a representative value for Case I, h* = 1.8 mm. For Case 
IV, where the normal stress is ~ 1 MPa and A-UsUp ranges from about 10 m/s to 100 m/s, 
the minimum value of h* is 8.1 mm and for much of the velocity range is significantly 
greater than that. It is worth noting that it is for Case IV that a supersonic pulse is 
obtained. For the other cases, h* typically has values between those for Cases I and IV. 
Several calculations were carried out to assess the mesh dependence of the results. For 
Case I (crack-like mode) and for Case III (pulse-like mode), calculations were carried out 
with the mesh spacing halved along the interface (0.15 mm in the uniform region) so that 
the number of elements was doubled. For Case I, with the finer mesh the transition region 
at the sliding tip consisted of four elements (sixteen integration points) as opposed to two 
elements (eight integration points) with the coarser mesh. The values of A-UgUp obtained 
from the two calculations nearly coincide. The mean values of the curves of Vnp versus 
time computed using the finer mesh also essentially coincide with the curves in Fig. EJ but 
the oscillations about the mean differ somewhat. In addition, the effect of mesh size was 
further explored for a crack-like mode and a pulse like mode using frictional parameters 
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that differ from those used in Cases I and II. The reference mesh shown in Fig. ^p and 
two other meshes were used for this. All the meshes have the smallest elements near 
the interface and the impact edge as for the reference mesh in Fig. ^p. The coarse mesh 
consisted of 1, 152 elements and 4, 812 degrees of freedom, with the smallest element size 
2.1 mm X 1.0 mm. The fine mesh consisted of 56,832 elements and 229,364 degrees of 
freedom, with the smallest element size 0.075 mm x 0.1 mm. The slip rate at a point along 
the interface was plotted as a function of time for both meshes. In the calculation with the 
crack-like sliding mode, the results from all three calculations are in very good agreement 
for the propagation speed except that increased fluctuations around a steady slip rate 
occur with the coarse mesh. For the pulse-like mode calculation, there is only a small 
difference between the calculations using the three meshes for the predicted magnitude 
and timing of the pulses. 

The values used for the interface stiffnesses Cn and Cs in eqs. © and ((Tj) play a signif- 
icant role both for symmetric loading (with no sliding) and for asymmetric loading (with 
frictional sliding). For symmetric loading. Fig. [71 the width of the head wave emanating 
from the loading wave front decreases with increasing interface stiffness. Because in the 
experiments the width of this head wave expands with decreasing applied load, the sym- 
metric loading geometry can be used to calibrate the stiffness values. The values used 
in most of the calculations, Cn = 0.3 MPa/m and Cg = 0.1 MPa/m, give rather good 
agreement for symmetrically loaded specimens with Eq in the range 5 — 40 MPa. 

To assess the dependence of the sliding mode on the interface stiffness values employed, 
asymmetric loading calculations were carried out using various values of C„ and Cs- For 
given values of Eq and l^mp, the frictional sliding mode that occurs depends on C„ and Cg- 
With Cn = 3.0 MPa/m and Cg = 1.0 MPa/m, only crack-like behavior was obtained for 
the range of values of Eg and V^mp used in Fig. UHl On the other hand, various calculations 
were carried out with C„ = 0.03 MPa/m and Cg = 0.01 MPa/m and a train of pulses was 
not obtained, although as shown for Case IV pulses coalescing to form a crack-like sliding 
mode did occur. Also, pulse steepening became more pronounced when Cn and Cg were 
decreased. 

The elasticity of the interface plays a fundamental role because the frictional slip 
A-Usiip is an internal variable determined by the frictional law (analogous to the plastic 
strain rate in plasticity theory). The observable quantity is the jump in displacement rate 
across the interface, Au^. The pulses occur in Aitshp not in Awg and, hence, according to 
our numerical results, are not directly observable in terms of a velocity jump distribution 
along an interface. What is observable is a pulse-type traction rate distribution fFigs.lTUb 
and 112b). At a given spatial location, the oscillatory nature of the pulse-like modes is 
seen in the time evolution of A-Ug (Fig. I^Uk). 

Calculations modeling dynamic slip in earthquakes, e.g. Cochard and Madariaga 
(1994, 1996), Perrin et al. (1995), Beeler and Tullis (1996) (see also Ben-Zion, 2001; 
Rice, 2001; Nielsen and Madariaga, 2003), exhibit both crack-like and pulse-like sliding 
modes. In analyzing conditions that set the mode of sliding, Zheng and Rice (1998) 
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considered the special case of identical material properties on both sides of the slip surface. 
Specifically, Zheng and Rice (1998) considered elastic half-spaces that are subject to a 
uniform shear stress Tq outside a perturbed region where slip nucleates. The criterion 
delineating between crack-like and pulse-like sliding involves the magnitude of Tq, the 
elastic moduli and shear wave speed of the half-spaces, the sliding speed and the sliding 
speed dependence of the interface friction law. The criterion of Zheng and Rice (1998) 
predicts that low values of Tq favor pulse-like sliding and higher values crack-like sliding. 

The circumstances considered here are quite different from those in Zheng and Rice 
(1998). In particular, it is important to note that the parameters in Fig.^lare the initial 
compressive stress and the impact velocity. In this parameter space, a clear demarcation 
is found between the regime where the crack-like and pulse-like modes occur; a low initial 
compressive stress and/or a low impact velocity favor the crack-like mode. 

The experiments and calculations exhibit a variety of common features including: 

(i) Direct calculation from the time history of the position of the sliding tip and Mach 
lines in the contours of cti — (J2 indicate that the sliding tip travels at a speed greater 
than ^/2cs and that a trailing pulse travels faster than q. 

(ii) The fringe patterns in front of the sliding region have similar characteristic shapes. 

(iii) Crack-like sliding occurs at low impact velocities and a sliding mode with a pulse is 
obtained at higher impact velocities. 

(iv) Two shock waves eventually merge into a single shock wave as the trailing shock 
catches up to the leading one. The tip of the trailing pulse travels at a speed 
exceeding c/. 

The calculations exhibit two sliding modes with two approaching shock waves. One 
involves trailing pulses traveling faster than the leading pulse that eventually coalesce to 
form one sliding region (Figs. IT31 and IT^ . The other involves a steepening and narrowing 
pulse followed by a crack-like sliding region (Figs. IT7I and ITHj) . 

There is some evidence in other contexts regarding the range of sliding modes that 
emerge from our calculations, more specifically regarding the pulse and pulse train modes 
of frictional sliding. Rubio and Galiano (1994) carried out experiments on sheared gels 
sliding along smooth glass surfaces and observed sliding via propagation of a quasi-periodic 
pattern of sliding zones of finite width separated by non- moving regions, having a propaga- 
tion speed of the order of q of the gel. This is consistent with our observations of relative 
speeds of pulse propagation that are three order of magnitude higher than the imposed 
impact velocities. The experiments of Rubio and Galiano (1994) motivated an analysis 
by Caroli (2000) of sliding of a viscoelastic solid on a rigid substrate with a rate and state 
friction law that showed a pulse train. Sliding between a rectangular Polyurethane slab 
and a compressed Araldite plate was studied experimentally by Mouwakeh et al. (1991). 
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The velocity measured at a fixed spatial point was strongly oscillatory. Also, the mea- 
sured traction exhibited oscillations. The traction versus time curve in Mouwakeh et al. 
(1991) resembles one of the oscillatory curve in Fig.l2Th. The velocity versus time curve is 
qualitatively similar to the intermediate curves in Fig. 121b. However it is unclear whether 
Ausiip or Ails is plotted in Mouwakeh et al. (1991, Fig. 5a). 

Fiber pull-out experiments exhibit a variety of complex frictional phenomena, e.g. 
Tsai and Kim (1996) and Li et al. (2002). In particular, Tsai and Kim (1996) observed 
three frictional sliding modes: steady sliding of the entire contact surface, stick-slip sliding 
of the entire contact surface, and sliding through the generation of concentrated sliding 
pulses. In one set of experiments, with sliding along the entire contact surface, the pull- 
out force was steady for a low interface pressure and was oscillatory for a high interface 
pressure, analogous to the behavior in our calculations where a crack-like sliding mode 
occurs for low values of Sq and a pulse-like mode for high values of Sq. However, Tsai 
and Kim (1996) found that the oscillation amplitude decreased with increasing pull-out 
speed. This trend, also found in Baumberger et al. (1994) and Povirk and Needleman 
(1993), under quasi-static loading conditions, is opposite to what occurs in our analyses 
under dynamic loading conditions. 

Brune et al. (1993) and Anooshehpoor and Brune (1994) observed two modes of 
sliding along the interface between identical foam rubber blocks; one involving wrinkle- 
like waves and the other a crack-like sliding mode. Characteristics of the wrinkle-like wave 
mode in these experiments are similar to those obtained for slip-pulses in our calculations; 
for example, an increasing normal component of particle motion with increasing normal 
stress which is analogous to the slip pulse amplitude increasing with Sq in our calculations. 
Interestingly, it was found that the frictional heating did not change with shear stress in 
the wrinkle-like wave mode whereas it increased linearly with shear stress in the crack-like 
mode (see also Andrews and Ben-Zion, 1997, Fig. 12). In experiments on sliding of two 
dissimilar blocks of foam rubber, Anooshehpoor and Brune (1999) observed that sliding 
occurred through a series of multiple opening pulses. They also found that the sharpness 
and the amplitude of the pulses either increased or decreased with time and was not steady. 
Such opening pulses, often referred to as Schallamach (1971) waves, have been observed 
experimentally in various bimaterial systems (e.g. Schallamach, 1971; Anooshehpoor and 
Brune, 1999) and in simulations (Coker et al., 2003). There is evidence for slip pulses, on 
which attention is focused here, in the experiments of Mouwakeh et al. (1991) and Tsai 
and Kim (1996). 

Figure El shows that the displacement distribution from a pulse train slip mode and 
from a crack-like slip mode are quite similar. This might be the reason that there is 
ample experimental evidence of slip pulses in soft materials such as gels, whereas there 
is no conclusive experimental evidence for slip pulses in engineering materials where the 
displacement magnitudes are very small compared to those in soft materials. However, 
the shear traction distribution along the interface in Fig. El does provide an indicator of 
the slip-mode. The variation of the shear resistance and the step increase in displacement 
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in Fig. E]is similar to the response seen in rock friction experiments (TuUis, 1996) and in 
fiber pull-out experiments (Tsai and Kim, 1996). Since the resulting slip from a train of 
pulses is hard to distinguish from that occurring with the crack-like slip mode, our results 
suggest the possibility that slip pulses may be more common than previously appreciated. 
It remains to be determined whether there is a significant difference between the heat 
generation associated with the pulse-train mode and that associated with the crack-like 
mode. 

6 Conclusions 

Frictional sliding between two Homalite plates subject to a compressive stress was investi- 
gated both experimentally and computationally. Sliding under shear loading was induced 
by impacting one plate parallel to the interface. The experiments yield information on 
the sliding speed and the stress state in the material using dynamic photoelasticity. A 
plane stress initial/boundary value problem is analyzed using a framework where sliding 
is modeled by a rate- and state-dependent frictional interface constitutive relation. 

• The experiments exhibit both pulse-like and crack-like modes of sliding. 

• In the experiments, the speed of the leading pulse or the leading edge of the crack- 
like sliding region ranges from somewhat above a/2cs to q. This propagation speed 
increases with increasing impact velocity and decreasing compressive stress. A speed 
exceeding q was seen for a trailing pulse. 

• A variety of frictional sliding modes were obtained in the calculations, depending 
on the initial compressive stress, the impact velocity and the interface characteriza- 
tion: a crack-like mode; a pulse-like mode with well-separated pulses that increase 
in amplitude; and a train of pulses that propagate with an essentially constant am- 
plitude. In addition, combinations of these modes occurred as well as transitions 
between modes. This variety of sliding modes emerges even though there is no 
elastic mismatch across the interface. 

• The slip resulting from the pulse-train mode and that resulting from the crack-like 
mode are hard to distinguish. 

• In all calculations the speed of the leading pulse or the leading edge of the crack- 
like sliding region exceeds v^c^ and is close to (or slightly exceeds) q. As in the 
experiments, trailing pulses with a speed exceeding q are found. 

• The elasticity of the interface plays a significant role in setting the mode of sliding. 

• The range of sliding modes obtained appear to be generic, arising in a wide variety 
of configurations and applications, and at a wide variety of size scales. 
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Appendix: Time integration of the friction law 

The equations are integrated in time using of a rate tangent method (Peirce et al., 1984) and 
automatic time step control. The shding velocity, AusUp, in the interval [t,t + dt], is expressed 
as a linear combination of its values at t and t + dt via 



Ausiip = (1 - 7)A-Usiip(t) + 7Ansiip(i + dt) 

where 7 is a chosen integration parameter that can range in value from to 1. 
A first order expansion of AttsUp in terms of T^ , 61,62, and ^o gives 
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The integral along the interface in eq. (^ is written as the sum over linear displacement 
interface elements. All field quantities and constitutive state variables are known at time t. 
Also, the values of Am„ and Aiig have been determined from solving eq. Q. Thus, An„ and 
Aus are know at time t + dt. Then, at a given integration point along the interface: 

1. Calculate Aiig = [Aus{t + dt) - Aus{t)]/dt and Aiin = [Au„(t + dt) - Aun{t)]/dt. 

2. If Aun{t + dt) > Tiq/Cji then set r„ = T^ = and go to the next integration point. 

3. Using the values for the state variables at time t, calculate g{6o)- 

4. Compute 

^^ 9{eo){Oi + e2)- ^^^^ 

5. If /3 < 1 then Aii^up = 0. Go to[7| 

6. If /3 > 1 calculate Cf " and Rp in eq. ^ . 

7. Update the tractions and state variables via 

r„fc+i = T^ + fndt (33) 

j.fc+1 ^ T^ + tsdt (34) 

^f+i = e^ + Oidt, i = 0,1,2 (35) 
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Property 


Value 


(J'S 


0.6 


l^d 


0.5 


Vo (m/s) 


100 


Vi (m/s) 


26 


B 


4.60 


P 


1.2 


m 


5 


Lq (m) 


20 X 10-6 


Li (m) 


20 X 10-6 


L2 (m) 


20 X 10-6 


C (MPa-i) 


0.7 


D (MPa-^) 


0.3 



Table 1: Friction parameter values used in the finite element calculations. 
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Figure 1: (a) Geometry and loading configuration used in the plane stress finite element 
calculations, (b) The finite element mesh used in the calculations (56, 580 degrees of 
freedom) . 



0.65 



0.55 




0.45 



0.65 
0.6 

I- 

0.5 

0.45 

0.4 



Au ,. (m/s) 
slip ^ ' 



10 m/s 


1 00 m/s 30 m/s 


- 


\ 




V r 


1 


1 1 1 1 



6 

(b) 



8 



10 

Au , /L 
slip 



12 



Figure 2: (a) Comparison of the steady-state coefficient of friction, figg, as a function of 
Awsiip for eq. (fT^ (dashed hne) and for eq. (jTHj) using the properties in Table [T] (sohd 
hne). (b) The effect of an abrupt change in A-UsUp on the apparent coefficient of friction 
Ts/Tn for the rate- and state-dependent friction relation used in the calculations. 
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Figure 3: Variation of the shear traction, Tg, (sohd hne) as a function of AugUp for a step 
jump in compressive normal traction, T„, (dashed hne) at constant frictional shding rate 
A-Usiip for the rate- and state-dependent friction relation used in the calculations. 
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Figure 4: Experimental isochromatic fringe patterns from a dynamic friction experiment 
on Homalite subject to a static compressive stress of 9.4 MPa and an impact velocity of 
32.7 m/s. In the inset lines are drawn to highlight Mach lines. The field of view is 130 
mm in diameter. 
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Figure 5: Experimental isochromatic fringe patterns from a dynamic friction experiment 
on Homalite subject to a static compressive stress of 9.4 MPa and impact velocity of 42 
m/s, at (a) t = 40/is, (b)t = 48/is, (c)t = 60fis. In the inset a line or lines are drawn to 
highlight Mach lines. The field of view is 130 mm in diameter. 
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Figure 6: Propagation speeds, K-upt, versus time for the five cases analyzed. The symbols 
are the experimentally measured propagation speeds from the Mach cone angle with 
Eo = 9.4 MPa and M^p = 42 m/s (Fig. \^. 
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Figure 7: Isochromatic fringe pattern due to symmetric loading with respect to the in- 
terface, (a) Numerical contours at t = 16 /is for Sq = 10 MPa, \4mp = 5 m/s, with the 
configuration analyzed shown in the inset, (b) Experimental isochromatic fringe patterns 
for So = MPa, Vimp = 58 m/s. In (b) the field of view is 130 mm in diameter. 
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Figure 8: Distributions along a portion of the interface for Case I - Sq = 6 MPa, V^mp = 2 
m/s. (a) The frictional shding rate, A-Usiip, and the shear traction, T^ at t = 22 /xs and at 
t = 32 /is. The symbols mark values at integration points along the interface to illustrate 
the resolution of the discretization, (b) The shear traction increment, dTg = Tsdt, and 
the velocity jump across the interface, A-Ug, at t = 32 /is. 
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Figure 9: Contours of cJi — (72 (isochromatic fringe patterns) at t = 32 /xs for Case I. 
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Figure 10: Distributions along a portion of the interface for Case II - So = 30 MPa, 
Vimp = 2 m/s. (a) The frictional sliding rate, A-UgUp, and the shear traction, T^ at 
t = 32 /is and at t = 34 /xs. The symbols mark values at integration points along 
the interface to illustrate the resolution of the discretization, (b) The shear traction 
increment, dT^ = Tgdt, and the velocity jump across the interface, Aiis, at t = 34 fis. 
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Figure 11: Contours of cxi — a^ (isochromatic fringe patterns) at t = 34 /xs for Case II. 
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Figure 12: Distributions along a portion of the interface for Case III - Sq = 10 MPa, 
^imp = 20 m/s. (a) The frictional shding rate, A-UsUp, and the shear traction, Tg at t = 32 
/is and at t = 33 /is. (b) The shear traction increment, dT^ = Tgdt, and the velocity jump 
across the interface, Aiig, at t = 32 /is. 




Figure 13: Contours of cJi — (72 (iso chromatic fringe patterns) at t = 32 /xs for Case III. 
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Figure 14: Distributions of the accumulated frictional sliding, A-UsUp = / '^UsUpdt, the 
displacement jump across the interface, Aus and the apparent coefficient of friction, /iapp = 
Ts/Tn, across the interface, (a) Case I - Sq = 6 MPa, Vimp = 2 m/s. (b) Case III - Sq = 10 
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Figure 15: Distribution of the frictional sliding rate, A-UsUp for Case IV - Sq = 
Mmp = 10 m/s. Also, for Case IV C„ = 0.03 MPa/m and C, = 0.01 MPa/m. 
/is. (b) t = 32 /is. (c) t = 48 /is. 
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Figure 16: Contours of o"i — o"2 (isochromatic fringe patterns) for Case IV. The contour 
lines are drawn to focus on details of the distribution for X2 > 0. (a) t = 22 fis. (b) t = 32 
/is. (c) t = 48 /is. The inset in (b) indicates the Mach line orientations. 
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Figure 17: Distribution of AiisUp at t = 32 /is, t = 38 fis and t = 44 fis and T, at t = 44 
/is along a portion of the interface for Case V - Sq = 40 MPa, Vimp = 2 m/s. 
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Figure 18: Contours of ai — a2 (isochromatic fringe patterns) for Case V. (a) t = 32 fis. 
(b) t = 38 fis. (c) t = 44.6 /xs. The dark lines indicate the location and width of the 
pulses in AusUp in Fig. [T3 
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Figure 19: Dependence of the mode of sliding on the initial compressive stress Sq and the 
impact velocity Vimp- Transition cases where mixed crack-like mode and pulse-like modes 
occur are shown as triangles. 
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Figure 20: Effect of varying impact velocity, Mmp, on tlie frictional sliding mode for Sq 
fixed at 10 MPa. (a) Aiig, (b) A-UgUp, (c) T^ and yUapp = Ts/Tn- 
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Figure 21: Effect of varying initial compressive stress, So, on the frictional sliding mode 
for \4inp fixed at 2 m/s. (a) Aiis, (b) Awsiip, (c) T^ and /iapp = Ts/T„. 



